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RADAR  ABSORBER  MEASUREMENT  TECHNIQUES 
AT  FREQUENCIES  ABOVE  20  GHz 

N.  S.  Nahman,  C.  M.  Allred,  J.  R.  Andrews 

C.  A.  Hoer  and  R.  A.  Lawton 


New  methods  for  implementing  automatic  permitti vity/permeabi 1 ity  measurements  of  radar 
absorber  materials  for  applications  primarily  above  20  GHz  and  into  the  millimeter  wave 
region  to  about  100  GHz  are  discussed.  A brief  review  of  the  state-of-the-art  of  dielectric 
and  magnetic  material  measurements  is  given.  Automated  time  domain  and  frequency  domain 
methods  are  considered  including  time  domain  automatic  network  analyzers,  correlation 
measurement  systems  using  noise  signal  sources,  and  six-port  network  analyzers. 

Key  words:  Dielectrics;  frequency  domain  measurements;  magnetic  materials; 

microwaves;  millimeter  waves;  permeability  measurements;  permittivity  measurements; 
radar  absorbers;  time  domain  measurements. 


1.  INTRODUCTION 
1.1  Review  and  Perspectives 

The  purpose  of  this  report  is  to  describe  new  methods  for  implementing  automatic  mu-epsilon  (p-e) 
measurements  of  radar  absorber  materials  for  applications  primarily  above  20  GHz  and  into  the 
millimeter  wave  region  to  about  100  GHz. 

An  adequate  view  of  the  state-of-the-art  for  complex  p and/or  e measurements  can  be  obtained  from 
the  technical  journals;  a typical  collection  of  technical  papers  is  given  in  references  [1.1-1.10], 
Because  of  the  need  to  rapidly  obtain  accurate  data  for  these  electromagnetic  parameters  over  broad 
ranges  of  frequency  and  in  extreme  and  changing  environments,  automated  measurement  methods  are 
requi red. 

Most  of  the  cited  present  day  techniques  are  narrowband  techniques  [1.1-1. 6]  with  the  exception  of 
Dispersive  Frequency  Domain  Spectroscopy  (DFDS)  [1.7, 1.8]  and  Time  Domain  Spectroscopy  (TDS) 

[1.9,1.10]  the  latter  two  techniques  employ  broadband  spectral  sources  in  the  form  of  incoherent  and 
pulsed  sources,  respectively.  DFDS  is  employed  in  the  optical /infrared  regions,  while  TDS  is  used  in 
the  electrical  regions  including  the  higher  microwave  frequencies. 

Regarding  terminology,  in  the  electrical  region  (wavelengths  longer  than  those  of  the  infrared 
region)  noise  sources  and  correlators  comprise  a system  equivalent  to  the  optical/infrared  DFDS  system. 
Both  systems  are  based  upon  correlation  functions  (observables)  and  their  Fourier  transforms. 

Similarly,  TDS  is  based  upon  time  domain  observables  and  their  Fourier  transforms.  In  this  case  the 
signal  sources  are  repetitive  and  of  finite  duration,  and  sometimes,  periodic. 

Recently,  new  single-frequency  measurement  techniques  based  upon  six-port  or  n-port  devices  have 
been  developed  for  precision  measurements;  these  techniques  are  amenable  to  calibration  using  fewer 
standards  but  require  on-line  computer  data  acquisition  and  processing.  However,  because  of  tire 
computer  tie-in,  the  frequency  can  be  rapidly  stepped  across  a frequency  band  to  provide  broadband 
frequency  domain  measurements.  Consequently,  such  measurement  methods  could  be  used  to  determine 
complex  p-e  data  across  a broad  frequency  band. 


1.2  Report  Content  and  Organization 


This  report  is  divided  into  six  chapters,  the  first  one  being  the  present  introduction.  Four 
chapters  are  devoted  to  time  domain  spectroscopy,  correlation  measurement  systems  with  noise  signal 
sources,  and  frequency  domain  automatic  systems  using  six-port  devices.  References  appear  at  the  end 
of  each  respective  chapter.  The  remaining  chapter  includes  a summary  and  recommendations  for  future 
research  and  development. 

Each  of  these  three  new  techniques  for  y-e  measurements  are  practical  techniques.  Each  one  has 
its  own  advantages  depending  upon  various  considerations.  In  establishing  the  feasibility  of  these  new 
methods,  the  MBS  contributors  have  developed  the  analyses  presented  herein  and  also  have  drawn  upon 
their  respective  experiences  in  research  and  development  on  rel ated  measurement  techniques  in  the  time 
domain,  correlation,  and  microwave  technical  areas.  Typical  factors  such  as  present  measurement 
capabilities,  cost,  staff  experience,  etc.,  will  most  probably  be  the  controlling  ones  for  the 
selection  of  one  method  over  another.  Also,  in  one  frequency  band  a given  technique  might  be  more 
attractive  than  another  for  the  very  same  reasons  or  combinations  of  the  typical  factors  cited  above. 

In  short,  this  report  does  not  say  which  method  should  be  used  but  does  present  inclusive 
recommendations  for  future  technical  work.  Also,  the  report  provides  a technical  basis  for  discussions 
between  the  U.S.  Air  Force  (UoAF)  Avionics  Laboratory  and  MBS. 

Chapter  2 is  concerned  with  the  time  domain  automatic  network  analyzer  (TDANA)  as  applied  to 
p-e  measurements  and  discusses  two  topics:  (1)  extension  of  present  baseband  pulse  techniques  to 

higher  microwave  frequencies,  and  (2)  pulsed  carrier  techniques  for  microwave/millimeter  wave 
applications.  Typical  time  domain  systems  are  the  NBS  Automatic  Pulse  Measurement  System  (APMS)  and 
the  USAF  Avionics  Laboratory  Time  Domain  y-e  Measurement  System. 

Because  the  sampling  process  and  the  sampling-head  structure  fundamentally  limit  the  temporal 
resolution  (and  correspondi ng  bandwidth)  of  sampling  oscilloscope  systems,  improved  sampling  heads 
could  provide  increased  temporal  resolution.  Chapter  3 discusses  the  improvement  of  sampling-head 
hardware  using  photoconductor  gating  structures.  This  concept  was  invented  at  MBS  under  partial 
support  of  the  USAF.  The  patent  rights  have  been  assigned  to  the  USAF  [1.11]. 

Chapter  4 presents  the  fundamentals  of  using  noise  signal  sources  for  characterizing  two-port 
networks.  The  desired  broadband  network  or  material  temporal  parameters  are  extracted  by  correlation 
measurements,  and  the  correspondi ng  frequency  domain  parameters  are  obtained  from  the  Fourier 
transforms  of  the  correlation  (temporal)  functions. 

Chapter  5 introduces  frequency  domain  automatic  systems,  using  six-port  devices,  which  are 
basically  single-frequency  techniques  but  are  highly  amenable  to  automation  and  calibration.  This,  in 
turn,  enables  a sequence  of  rapid  measurements  to  be  made  over  a frequency  band. 

Finally,  Chapter  6 contains  a summary  of  the  body  of  the  report's  discussion  of  the  major  factors 
involved  in  implementation  of  new  methods,  and  recommendations  for  future  technical  work. 
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2.  THE  TIME  DOMAIN  AUTOMATIC  NETWORK  ANALYZER 


James  R.  Andrews 
2.1  Introduction 

Time  domain  techniques  have  been  applied  quite  successfully  to  the  measurement  of  frequency  domain 
parameters  at  frequencies  up  to  X-band  (12  GHz)  [2.1-2.10].  These  include  the  complex  S-parameters  for 
electrical  networks  and  material  parameters  [2.4,2.11,2.12]  such  as  the  complex  permeability  and 
permittivity  (y  and  e)  of  dielectrics  and  ferrites.  The  recent  review  paper  by  Andrews  summarizes  the 
techniques  and  state-of-the-art  [2.13].  The  common  features  of  the  time  domain  automatic  network 
analyzer  (TDANA)  (fig.  2.1)  as  used  by  the  various  researchers  include:  (1)  a very  fast  pulse 

generator,  (2)  a wideband  sampling  oscilloscope,  and  (3)  a minicomputer.  The  USAF  Avionics  Laboratory 
y-e  measurement  system  is  an  example  of  such  a system  [2.14-2,16]. 

The  pulse  generators  used  have  typically  been  either  tunnel  diode  step  function  generators  or  step 
recovery  diode  (SRD)  impulse  generators.  A step  recovery  diode  generator  typically  produces  a 5 V, 

60  ps  impulse  with  a useful  spectrum  up  to  8 GHz.  A tunnel  diode  typically  produces  a 1/4  V step  with 
a 20  ps  transition  time  and  a useful  spectrum  up  to  12  GHz.  Recently,  measurement  capabilities  have 
been  extended  to  the  Ku-band  (12.4-18  GHz)  using  improved  pulse  generators  which  generate  a short  rf 
burst  [2.17]. 

The  state-of-the-art  in  sampling  oscilloscopes  is  "A"  which  has  a 20  ps  transition  time  and  a dc 
to  18  GHz,  3 dB  bandwidth.  L Another  sampling  oscilloscope  often  used  is  the  "B".  It  has  a 30  ps 
transition  time  and  a dc  to  11.5  GHz,  3 dB  bandwidth. 

A minicomputer  serves  many  functions  in  the  operation  of  a TDANA,  including:  (1)  measurement 

system  controller,  (2)  data  acquisition,  (3)  data  storage,  (4)  digital  averaging  to  improve  signal-to- 
noise  ratio,  (5)  Fourier  transform  computation  to  convert  the  time  data  to  the  frequency  domain,  (6)  com- 
putation of  the  scattering  parameters  S..  or  other  desired  parameters,  and  (7)  presentation  of  the  results 
in  a useful  tabular  or  graphical  form  to"1  the  operator  via  a cathode  ray  tube  (CRT)  terminal. 

When  one  is  faced  with  making  S-parameter  or  y-e  measurements  in  the  millimeter  wave  region,  one 
question  which  arises  is  whether  the  TDANA  technique  can  be  extended  to  these  frequencies.  The  answer 
is  a qualified  "yes."  This  chapter  will  deal  with  this  question.  The  basic  TDANA  arrangement  with  a 
shorter  rf  pulse  burst  and  the  present  HP  sampler  can  probably  be  extended  to  the  K-band  (18-26.5  GHz) 
with  some  difficulty;  the  next  two  sections,  2.2  and  2.3,  will  address  this  topic.  Above  K-band, 
alternate  techniques  will  no  doubt  have  to  be  used.  At  the  present  time,  the  most  promising  appears  to 
be  that  of  heterodyning  a baseband  pulse  up  to  the  millimeter  band  of  interest  and  then  synchronously 
heterodyning  a response  signal  back  to  baseband  for  subsequent  measurement;  section  2.4  will  discuss 
this  topic. 


-'-In  order  to  avoid  influencing  the  competitive  positions  of  manufacturers  of  commercial  equipment 
we  have  not  identified  these  products. 
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2.2  Extension  of  Present  TDANA  Techniques 


The  first  question  to  be  faced  is,  "Will  the  "A"  sampler  work  at  millimeter  frequencies?”  The 
manufacturer  specifies  the  transition  time  to  be  « 20  ps  and  the  bandwidth  as  dc  to  greater  than  18  GHz. 
"A-B"  and  "A-C"  samplers  owned  by  NBS  do,  in  fact,  have  transition  times  of  a 20  ps  when  they  are  properly 
adjusted.  An  NBS  Guest  Worker,  S.  M.  Riad,  as  part  of  his  Ph.D.  dissertation,  mathematically  modeled  the 
older  "A-A"  sampler  [18].  The  manufacturer  specified  the  bandwidth  of  the  "A-A"  to  be  dc  to  > 12.4  GHz 
and  the  transition  time  to  be%28  ps.  Riad's  model  predicted  the  step  response  transition  time  to  be 
27.5  ps.  Figure  2.2  shows  Riad's  computed  impulse  response  for  the  "A-A".  The  impulse  response  duration 
(50%)  is  24.1  ps.  Figure  2.3  shows  the  frequency  response  computed  from  Riad's  impulse  response.  While 
it  is  not  mathematically  rigorous,  an  approximation  to  the  "A-B/C"  response  can  be  made  by  scaling 
figure  2.2  by  20/28  and  figure  2.3  by  18/12.4.  From  this  extrapolation,  the  "A-B/C"  amplitude  response 
may  be  estimated  as  follows: 


Table  I.  "A-B/C"  Predicted  Frequency  Response 
Ampl  i tilde 

Response  Frequency 


- 3 

dB 

18 

GHz 

- 6 

dB 

15 

GHz 

- 10 

dB 

35 

GHz 

- 14 

dB 

42 

GHz 

- 20 

dB 

50 

GHz 

Figure  2.4  is  offered  as  further  evidence  that  the  "A-B/C"  will  function  at  millimeter  frequencies. 

This  figure  shows  an  actual  45  GHz  sine  wave  that  was  observed  on  an  NBS  "A-B"  during  the 
development  of  the  SHF  impulse  generator  built  for  the  USAF  Avionics  Laboratory  in  1977  [2.17]. 

Figures  2. 2-2. 4 demonstrate  that  the  "A-B/C"  sampler  will  operate  above  18  GHz.  Because  operation 
is  in  the  high  frequency  roll-off  region,  measured  waveforms  will  be  considerably  distorted.  However, 
when  making  ratio  measurements,  such  as  S21  or  Sn,  the  frequency  response  effect  of  the 
sampler  transfer  function  Hs(tu),  will  cancel  out.  In  theory,  Hs(co)  is  immaterial  in  ratio 
measurements.  In  practice,  it  has  a significant  effect.  If  Hs(cu)  is  down,  then  the  signal-to-noise 
ratio  in  the  measurement  is  reducede  and,  as  a result,  the  measurement  errors  are  increased. 

Although  the  "A-B/C"  sampler  response  is  seen  to  be  useful  to  at  least  50  GHz,  there  are  other 
problems  that  present  major  obstacles  to  measurements  above  18  GHz.  They  are  primarily  timing 
instabilities,  such  as  jitter  and  long-term  drift.  The  source  of  these  instabilities  is  in  both  the 
sampling  oscilloscope  time  base  and  the  pulse  generator.  Cans  [2.8]  and  Elliott  [2.19]  have  shown  that 
the  effect  of  signal  averaging  on  a waveform  that  has  time  jitter  is  equivalent  to  using  a low  pass 
filter.  For  stationary,  Gaussian  time  jitter  of  standard  deviation,  oj,  the  filter  is  of  the  form 


-1/2  (27Tfo.)2 

H (f ) a e 


(2.1) 


Table  II  lists  the  jitter  observed  on  several  different  NBS  sampling  oscilloscopes  and  pulse 
generators.  These  are  values  representati ve  of  the  current  state-of-the-art.  These  reported  jitter 
values,  t • are  human  observer  measurements  of  the  nominal  width  of  the  jitter  band  excluding  some 
outlying  samples.  They  are  related  to  crj  by 


a.  - 1/2  T..  (2.2) 

J J 


The  3 ps  value  of  oj  reported  by  Cans  [2.8]  was  obtained  from  a computer  controlled  measurement  of 
the  complete  jitter  probability  distribution.  Also  included  in  Table  II  is  the  equivalent  filter 
cutoff  frequency  (-3  dB)  due  to  the  jitter. 
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Table  1 1.  Jitter  Characteri sties 


Sampling  Oscilloscope 

Pulse  Generator 

Tj(ps) 

fC0(GHz 

"A-B"  sampler  with 
"C"  time  base  and 
"D"  vertical  amp. 

"E"  tunnel  diode  and 
NBS  bias  supply 

Tj  ^ 5-6  ps 
(human  observation) 
oj  = 3 ps  [2.8] 
(computer  meas.) 

44 

"A-B"  sampler  and 
"F"  sampling  plug- 
in 

"E"  tunnel  diode  and 
NBS  bias  supply 

v2  ps  (low  freq.  ripple) 
M ps  (Gaussian  jitter) 
v5  ps  (overall ) 

53 

"A-B"  and  "F" 

"E"  tunnel  diode  and 
NBS  improved  automatic 
TD  bias  supply 

^2  ps  (low  freq.  ripple) 
M ps  (Gaussian  jitter) 
^5-6  ps  (overall) 

48 

"A-B"  and  "F" 

"E"  tunnel  diode  and 
!'G”  bias  supply 

%2  ps  (low  freq.  ripple) 
^5  ps  (Gaussian  jitter) 
^6-7  ps  (overall) 

41 

"A-B"  and  "F" 

NBS  SHF  Impulse 
Generator  [2.17] 

v3  ps  (low  freq.  ripple) 
M ps  (Gaussian  jitter) 
^6  ps  (overal 1 ) 

44 

"B"  sampler  with 
"H"  time  base  and 
"I”  vertical  amp. 

"E"  tunnel  diode  and 
NBS  improved  automatic 
TD  bias  supply 

M ps  (Gaussian) 

66 

"B",  "H"  and  "I" 

NBS  SHF  Impulse 

v4  ps  (Gaussian) 

66 

Generator 


The  other  major  timing  error  is  due  to  long-term  drift  (ms  or  longer).  Drifting  is  potentially 
more  serious  than  jitter.  Drifting  has  many  different  causes,  including  minute  temperature 
fluctuations,  line  voltage  changes,  etc.  MBS  has  found  it  necessary  to  operate  its  scopes  and 
generators  continuously,  to  minimize  these  effects.  It  is  possible  to  devise  special  hardware  and/or 
software  techniques  to  minimize  these  long-term  drifts.  B.  Elliott  devised  a hardware  technique 
that  allowed  him  to  reduce  by  a factor  of  500  the  drift  of  a tunnel  diode  generator/”A-A"  sampler 
combination  [2.20].  The  overall  drift  stability  of  his  system  was  about  20  fs  (0.02  ps). 

Nicholson  developed  a hardware/software,  3-point  scanning  technique  to  detect  and  correct  drift 
errors  [2.21].  His  techniques  are  used  in  the  present  p-e  measurement  system  at  the  USAF  Avionics 
Laboratory  [2.14].  Both  Elliott's  and  Nicholson's  techniques  require  the  use  of  an  additional 
reference  pulse  directly  related  to  the  signal  to  be  measured.  This  reference  pulse  is  required  to  be 
positioned  at  a precise  location  in  the  time  window  of  observation.  These  techniques  are  useful  in  a 
TDANA  system  that  is  hardwired  into  a fixed  configuration  and  that  is  only  used  for  a single  type  of 
measurement.  To  date,  these  techniques  have  not  been  incorporated  into  the  NBS  TDANA  because  it  is 
used  for  a wide  variety  of  measurements  employing  many  different  pulse  generators  and  time  windows  from 
100  ps  to  several  ms. 

Another  effect  that  often  causes  timing  problems  is  the  60  Hz  ground  currents  between  the  various 
pieces  of  equipment.  For  ps  resolution,  even  minute  ground  currents  can  cause  noticeable  degradation. 

To  extend  the  direct  measurement  performance  of  the  NBS  oscilloscope  above  18  GHz  will  require 
several  modifications  including:  (1)  "F-A"  time  base  changes  to  reduce  jitter  and  60  Hz 

i nterference,  (2)  system  ground  loop  isolation,  and  (3)  incorporation  of  hardware/software  drift 
stabilization. 
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2.3  Millimeter  Wave  Pulse  Generator  Development 


As  noted  earlier,  the  present  baseband  pulse  generators  used  in  TDANA's  produce  a useful  spectrum 
up  to  18  GHz.  To  extend  beyond  18  GHz  will  probabably  require  the  use  of  the  same  technique  used  to 
cover  the  Ku-band,  namely,  the  generation  of  a short  rf  burst.  Figure  2.5  shows  the  waveform  and 
spectrum  of  an  ideal  rf  burst.  The  equation  for  the  waveform  is 

v (t ) = VQ[u(t  + t /2)  - u ( t - t /2)]  sin(2irf  t).  (2.3) 


The  spectrum  amplitude,  S(f),  is  given  by 


S(f)  = V t 
o o 


sin(irAft  ) sin[ir(2f  + Af)t  ] 
o c o_ 

(irAft  ) [n(2f  + Af)  t ] 

o co 


where 


Af  = f - f . 

c 


(2.4) 


(2.5) 


Several  techniques  can  be  used  to  generate  an  rf  burst.  NBS  Tech.  Note  699  discusses  in  detail 
several  of  these  techniques.  One  technique  is  to  gate  the  output  of  a cw  sine  wave  generator.  PIN 
diode  modulators  or  double  balanced  mixers  can  be  used  as  the  gates.  2 ns  or  slower  transition  times 
are  typical  for  PIN  diodes  while  1 ns  are  typical  for  double  balanced  mixers.  Another  technique  that 
NBS  has  found  useful  for  generating  high  power,  broadband,  microwave  rf  bursts  is  to  shock-exci te  a 
high  power,  traveling-wave  tube  amplifier  (TWTA)  with  an  extremely  narrow  (<  100  ps)  baseband  impulse. 
Figures  2.6  and  2.7  show  the  waveforms  obtained  from  C-,  X-,  and  Ku-band  TWTA's.  The  TWTA  acts  as  both 
a wide  bandpass  filter  (typically  octave  bandwidths,  2 to  4 GHz,  etc.)  and  a high  gain  amplifier  (30  dB 
is  typical).  The  output  is  a short  rf  burst  of  a few  rf  cycles  with  a frequency  near  the  low  end  of 
the  TWTA  passband.  Peak  power  comparable  to  the  output  power  rating  of  the  TWTA  can  be  obtained.  As 
seen  in  figure  2.7,  broadband  spectrum  amplitudes  of  the  order  of  80  dBpV/MHz  can  easily  be  obtained. 
The  C-band  TWTA  output  is  useful  (20  dB  range)  from  2 to  13  GHz.  The  X-band  TWTA  useful  range  extends 
from  4 to  15.5  GHz  with  the  exception  of  a sharp  notch  at  9.4  GHz.  The  Ku-band  TWTA  appears  to  be 
useful  from  approximately  10  GHz  to  16  GHz. 

A second  technique  for  generating  a microwave  rf  burst  is  to  mount  a step  recovery  diode  within  a 
piece  of  waveguide  [2.17,2.22].  Figures  2.8  and  2.9  show  the  waveform  and  spectrum  produced  by  the  SHF 
impulse  generator  recently  built  by  NBS  for  the  USAF  Avionics  Laboratory.  The  peak-to-peak  amplitude 
is  0.8  V.  The  center  frequency  is  12.5  GHz.  The  1/e  duration  is  1/4  ns.  The  spectrum  amplitude  is 
32  dByV/MHz,  + 7.5  dB , from  6 to  18  GHz.  The  step  recovery  diodes  used  in  the  SHF  impulse  generator 
were  the  fastest  (60  ps)  diodes  commercially  available.  They  could  also  be  used  in  K-band  waveguide 
(RG-53/U)  but  with  reduced  output  amplitude.  If  necessary,  the  amplitude  of  this  rf  burst  could  be 
increased  considerably  using  a high  power  TWTA. 

Techniques  other  than  the  step-recovery  (snap  -off)  transition  may  be  more  useful  at  K-band  and 
higher.  Other  examples  are  the  various  types  of  transferred  electron  devices  such  as:  avalanche, 

impatt,  trapett,  Gunn  diodes,  etc.  Figure  2.4  shown  earlier  is  an  example.  During  the  development  of 
the  NBS  SHF  impulse  generator,  the  driving  pulse  to  the  waveguide-mounted,  step-recovery  diode  was 
i nadvertantly  adjusted  to  exceed  the  reverse  breakdown  voltage  of  the  diode.  When  the  diode  started 
conducting  in  the  reverse  direction,  the  45  GHz  oscillations  of  figure  2.4  were  observed.  Various 
solid-state  oscillators  are  now  available  at  frequencies  up  to  at  least  200  GHz  [2.23].  Extremely  fast 
(<  100  ps)  on-off  modulation  of  some  of  these  oscillators  should  be  possible  with  some  modifications. 
Another  idea  that  should  be  investigated  is  the  use  of  a short  laser  pulse  to  trigger  on  an  oscillator 
diode. 
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2.4  Heterodyne  Techniques 


Heterodyne  techniques  have  been  used  since  the  earliest  days  of  wireless.  They  are  a 
straightforward  means  of  translating  from  one  frequency  to  another.  Heterodyning  will  permit  the 
extension  of  baseband  pulse  measurement  techniques  to  millimeter  wavelengths.  The  baseband  spectrum  of 
a test  pulse  generator  can  be  heterodyned  up  to  a high  frequency  and  passed  through  an  unknown  network 
or  material.  The  response  signal  can  then  be  downconverted  back  to  baseband  for  subsequent  measurement 
and  analysis.  This  section  will  discuss  the  advantages  and  disadvantages  of  various  modulation 
techniques  including:  amplitude  modulation  (AM),  frequency  modulation  (FM),  phase  modulation  (PM),  and 
single  sideband  (SSB).  Detection  techniques  discussed  include:  amplitude  detectors,  both  linear  and 

square  law,  frequency  discriminator,  phase  detector,  coherent  detector,  synchronous  detector, 
phase-locked  detector,  and  image  rejection  mixers. 

In  general,  a modulated  signal  may  be  described  as 


v(t)  = v (t)  sin[co(t)  t + (4(t)]. 

3 


(2.6) 


In  other  words,  we  can  modulate  either  the  amplitude  [ va (t ) ] , the  frequency  [a)(t ) ] , or  the  phase 
[«S(t)]  or  any  combination  of  the  above.  The  modulating  signal  vm(t)  can  be  of  any  form  such  as 
periodic,  aperiodic,  single  transient,  or  random.  For  the  purposes  of  time  domain  testing,  we  will 
hereafter  limit  the  discussion  to  periodic  signals.  The  periodic  signal  may  be  of  any  waveshape, 
including:  sine  wave,  impulse,  doublet,  rectangular  pulse,  trapezoidal  pulse,  rf  burst,  etc.  Being 

periodic,  it  can  be  completely  described  by  a Fourier  series  expansion  of  an  infinite  sum  of  sines  and 
cosines. 


v (t)  = ~ + 2 [A  cos(2nnf  t)  + B sin(27rnf  t)] 
m 2 , n on  o 

n=l 


(2.7) 


where  f0  is  the  fundamental  frequency  of  the  periodic  waveform.  The  frequency  coefficients  are 
gi ven  by 


A = r /v  (t)  cos  (2irnf  t)dt 
n T m o 


(2.8) 


2 T 

B =—  /v  (t)  sin(27rnf  t)dt 
n T m o 


(2.9) 


where  T = l/f0  is  the  period.  The  Fourier  series  can  also  be  written  in  complex  exponential  form. 


+«>  +j2irnf  t 

v (t)  = E D (nf  ) e 
m no 

n=-co 


(2.10) 


where 


dt  . 


(2.11) 


D 

n 


1 

T 


T -j27rnf  t 

/ v (t)e 
m 
o 


An  and  Bn  are  real  constants  while  Dn  is  a complex  constant.  Knowing  full  well  that  an 
arbitrary,  periodic  test  signal  can  be  represented  by  a series  of  sines  and  cosines,  we  will  hereafter 
limit  the  discussion  to  a single  term  in  the  Fourier  series.  Additional  terms  will  be  added  in  a real 
situation  but  would  only  serve  to  confuse  the  following  discussion. 

2.4.1  Amplitude  Modulation  - Double  Sideband 

Amplitude  modulation  (AM)  can  be  represented  as  shown  in  figure  2.10.  This  type  of  modulation  may 
be  described  by  the  equations: 


v (t)  = k.  v (t)  V cos  (co  t)  (2.12) 

rf  1 m c c 

v (t)  = V [1  + m cos  (co  t)]  (2.13) 

mm  m 


where  m is  the  modulation  index  and  k-j  is  the  ith  mixer  constant  (i  = 1 for  the  above). 


v ,(t) 

rf 


k.  V V 
1 o c 


= cos 


(“c*> 


cos 


(lilt)  + — COS  (co  t) 
l 2 u 


(2.14) 


a)  = co  - to  and  co  = co  + co 
l c m u c m 


This  leads  to  the  classical  result  with  AM  that  the  modulation  of  a single  carrier  frequency,  coc, 
by  a single  modulation  frequency,  (%,  creates  two  sideband  frequencies  spaced  co^  above  and  below 
the  carrier  frequency.  The  carrier  frequency,  coc,  term  can  be  considered  the  same  as  the  dc  term. 

For  further  Fourier  analysis,  it  can  be  dropped.  Actually,  if  a double  balanced  mixer  is  used,  the  dc 
component  of  the  modulation  is  removed,  and  m is  set  to  1,  then  the  carrier  will  be  suppressed  and  the 
output  will  only  consist  of  two  sidebands  (DSB).  Thus  eq  (2.14)  becomes: 

v -c(t)  = k V cos(co  t)  V cos(co  t)  (2.15) 

rf  1 m me  c 


V V 
m c 

vrf(t)  = k — — [cos  (a)  t)  + cos  (co^t)] 


(2.16) 


If  the  rf  test  signal  is  then  passed  through  an  unknown  network,  H(co),  the  output  is  given  by 


V (co)  = H(co)  V Jio). 
o rf 


(2.17) 


In  general , the  network  transfer  function  is  frequency  dependent  and  is  not  symmetrical  about  coc. 
Although  we  are  concerned  with  measuring  networks  or  materials,  this  same  problem  arises  in  short-wave 
radio  communication.  There,  the  ionosphere  creates  selective  fading  which  can  cause  severe  distortion 
of  the  received  signal.  In  general,  the  relative  amplitudes  and  phases  of  the  two  sidebands  will  be 
altered,  but  not  in  the  same  way  for  each  sideband.  Thus,  the  output  equation  must  reflect  this 
effect. 
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(2.18) 


vo(t) 

After  some  mathematical 

v (t) 
o 


kl  Vm  Vc 

r - [H  cos  (go  t + fa „)  + H cos  (u)  t + fa „)]. 
2 a i a u u a 


manipulation,  this  can  be  expressed  in  the  form 


k . V V joj  t ju)  t 

1 m c ..  „ u5,m  m c 

H Re  [ e e e 

2 ulm 


(2.19) 


where 


and 


u5,m 


[(H  cos 
u 


H cos 
l 


1/2 


+ (H  sin  fa  + H sin  fa  )] 
u u £ l 


(2.20) 


^u5,m 


tan 


H s i n 0 + H s i n 0 
u u H 5, 


H cos  fa  - H cos  fa 
u u i l 


(2.21) 


An  important  conclusion  to  note  from  the  above  is  that  in  general  we  have  two  measurable  quantities, 
namely  HU2,m  and  ^ ug,m , but  they  are  functions  of  four  unknown  variables,  namely,  Hu,  fau, 

H , and  fa^. 

The  next  question  that  needs  to  be  addressed  is  the  detection  of  the  rf  signal.  The  simplest 
technique  is  the  diode  envelope  detector  (fig.  2.11).  The  diode  detector's  output  is  proportional  only 
to  the  envelope  amplitude  of  the  rf  signal. 


vd(t)  01 


(2.22) 


A nonlinear  element,  typically  a semiconductor  diode,  is  used  to  perform  the  detection  process. 

For  large  signal  levels,  the  diode  output  is  linear  (p  = 1).  For  low  level  signals,  the  detected 
output  typically  follows  the  square  law  (p  = 2).  However,  in  practice,  neither  the  truly  linear  nor 
the  square  law  detector  is  really  achieved.  With  a proper  choice  of  input  level  and  circuit  constants, 
the  ideal  conditions  can  be  approached  closely. 

The  AM  (DSB)  signal  involves  four  unknowns,  and  the  simple  diode  detector  gives  only  one  measurable 
parameter.  Thus  the  use  of  AM  (or  DSB)  and  a diode  detector  is  not  the  proper  choice  for  rf  hetero- 
dyned network  measurements. 

The  synchronous  detector  (fig.  2.12)  is  a detection  scheme  that  permits  the  measurement  of  both 
the  amplitude  and  phase  of  an  rf  signal.  Examination  of  figure  2.12  shows  that  a synchronous  detector 
is  simply  a heterodyne  mixer  fed  with  a local  oscillator  identical  in  frequency  to  the  incoming  rf 
signal.  For  the  situation  shown  in  figure  2.12,  the  detector  output,  vscj,  is 


* ,(t)  = k v (t)  V cos  (go  t) 
sd  2 o 5,0  c 


(2.23) 


When  eq  (2.18)  is  used  for  v0(t),  vsc|  is  found  to  contain  terms  at  frequencies  (%  and 
2 goc  + GOffl • The  low  pass  filter  rejects  the  2ioc  + terms  with  the  result  that  vscj  is 
gi ven  by 
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(2.24) 


v ,(t) 
sd 


k.  k„V  V V 

1 2 m c 2.0  r..  , . . , . , „ 

- [H  cos  (a)  t + 0 ) + H cos  (a)  t - 0 )] 

4 u m u H m H 


Thus,  the  amplitude  and  phase  information  of  have  been  directly  translated  from  rf  frequencies 
to  baseband  frequencies.  From  eq  (2.24)  we  can  write  two  equations  (amplitude  and  phase  of  the 
modulation  frequency,  (%),  but  they  still  involve  four  unknowns,  H^,  0^,  Hu,  and  j5u. 

What  is  really  needed  is  some  technique  for  selectively  removing  one  of  the  sidebands  and  then 
measuring  the  network  response  to  the  remaining  sideband.  This  could  be  done  by  brute  force  filtering 
of  the  generated  DSB  signal.  However,  the  requirements  imposed  on  the  filter  would  be  quite  stringent. 
Also,  new  filters  would  be  required  each -time  one  wished  to  move  the  carrier  frequency. 

Another  technique  that  can  be  used  to  separate  the  information  from  the  two  sidebands  is  to  deal 
with  the  quadrature  nature  of  a DSB  signal.  Consider  again  the  DSB  case  described  earlier  in 
eqs  (2.15-2.18).  Equation  (2.18)  can  be  rewritten  in  the  form 

v (t)  = v.(t)  cos(u)  t)  - v (t)  sin  (cot).  (2.25) 

oi  c q c 


v-j(t)  and  Vq(t)  are  the  "in-phase"  and  "quadrature"  components,  respectively,  of  v0(t).  They 
are  given  by 


k V V 

v.  (t)  = — [H  cos(co  t + 0 ) + H cos  (oo  t - 0 )]  (2.26) 

i 2 u m u it  m H 

k V V 

v (t)  = - — C"  [H  sin(u)  t + 0 ) - H sin(co  t - 0 )]  . (2.27) 

q 2 u m u H m H 


v-j  and  Vq  can  be  detected  independently  if  two  synchronous  detectors  are  used  as  shown  in 
figure  2.13.  Note  that  the  LO  phase  for  the  quadrature  detector  is  shifted  by  90°.  The  in-phase 
detector  is  the  same  as  shown  in  figure  2.12  and  described  by  eqs  (2.23)  and  (2.24).  v -j ^ ( t ) is 
vscj(t)  as  given  by  eq  (2.24) . 


1/2 


k„V  v. (t) . 
2 Ho  i 


(2.28) 


The  quadrature  detector  output  is  given  by 

v (t)  = v (t)  k V [-sin(u)  t)  ]. 
qd  o 3 Ho  c 


(2.29) 


After  the  2wc  terms  are  filtered  out,  this  reduces  to 


Vt!  ' 1/2  k3 


Ho 


v (t). 
q 


(2.30) 


We  now  have  two  waveforms  from  which  we  can  measure  two  separate  magnitudes  and  phases.  All  that 
remains  is  to  find  a mathematical  separation  to  extract  the  four  unknowns  (H^,  0 Hu,  and 
«SU)  from  Vi<j  and  vqc|(t). 

The  Hilbert  transform  [2.24,2.25]  provides  the  key  to  unlocking  the  puzzle.  The  Hilbert  transform 
x(t)  of  a signal  x ( t ) is  given  by  the  principal  value  of  the  following  equation: 
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(2.31) 


i(t)  =H  Cx(t)]  = 1 f 7^  dT. 

IT  T 

_ oo 


x(t)  is  the  convolution  of  x(t)  with  1/t.  The  Hilbert  transform  can  be  thought  of  as  a network 
that  phase-shifts  the  positive  frequencies  by  -90°  and  the  negative,  by  90°.  The  Fourier  transform 
of  x(t)  is  X(f ) and  is  related  to  x(t)  and  X ( f ) by 


X ( F ) = -j(sgn  f)  X(f). 


(2.32) 


For  example,  if 


x (t ) = cos(u)t),  then  x(t)  = sin(ajt) 


(2.33) 


or  if 


t(t)  = sin(uit),  then  x(t)  = -cos(u)t) 


(2.34) 


Applying  these  results  to  the  detected  i and  q signals,  we  can  write  Vq<j(t)  as 


v ,(t)  = H[v  (t)]  = t k k V V V [-H  cos(<u  t + <b  ) + H cos(w  t - <t>  )].  (2.35) 

qd  qd  4 1 3 m c lo  u m u 2 m 2 


If  the  two  detector  mixers  are  identical 


k 


2 


= k 


(2.36) 


then  the  sums  and  differences  v-j<j  7 Vqj  will  give  the  separated  lower  and  upper  sideband 
components . 


Vl)  = ’to111 


qd 


( 7 k k V V V ) 
2 1 d m c to 


H cos  (co  t - 0 ) 
2 m2. 


(2.37) 


V,(t>  ' 


v ,(t)  = ( - k.kv  V V ) 
qd  2 1 d m c 20 


cos(w  t + <b  ) 
m u 


(2.38) 


Equations  (2.37)  and  (2.38)  were  the  sought  for  results.  They  show  that  we  can  use  DSB  if  we  use 
an  I and  Q synchronous  detector  and  further  process  the  quadrature  component  using  the  Hilbert 
transform,  v-j^t)  and  VqCj(t)  can  be  measured  directly  by  a TDANA  sampling  oscilloscope.  Then, 
the  Hilbert  transform,  VqC|(t),  and  the  sums  and  differences  v-jd(t)  ± Vqd(t)  can  be 
performed  by  the  TDANA  minicomputer.  It  is  also  possible  to  perform  these  operations  in  an  analog 
fashion  using  quadrature  hybrids  and  power  dividers/combiners.  An  I '5rT3'"Q_synchronous  detector  with 
these  additional  components  is  called  an  image  rejection  mixer  [2.26]. 


The  envelope  function,  ved,  can  also  be  expressed  in  terms  of  the  I & Q components. 


ved(t)  = (vid' 


, 2.1/2 

+ vqd  } 


(2.39) 


Figure  2.14  shows  a proposed  heterodyne  TDANA.  It  uses  DSB  signal  generation  and 
and  Q detection.  The  minicomputer  is  used  to  perform  the  necessary  Hilbert  transform. 


synchronous 
The  major 


I 


-12- 


bandwidth  limitation  in  this  scheme  is  the  available  i-f  bandwidth  of  the  double  balanced  microwave 

mixers.  DC  to  4 GHz  appears  to  be  the  state-of-the-art  at  present  [2.26].  With  this  setup,  a single 

TDANA  measurement  will  yield  data  over  an  8 GHz  frequency  span  (fc  + 4 GHz).  If  a 10  ns  time  window 

is  used,  then  the  frequency  spacing  of  the  resultant  data  would  be  100  MHz. 

2.4,2  Amplitude  Modulation  - Single  Sideband 

As  well  as  using  DS8 , one  should  also  consider  using  single  sideband  (SSB).  The  brute  force 
filtering  technique  of  SSB  generation  was  mentioned  earlier  and  dismissed  as  impractical  for  this 
application.  SSB  can  also  be  generated  using  the  phasing  technique  [2.24-2.29].  It  is  more  adaptable 
to  the  requirement  of  moving  the  carrier  frequency.  Figure  2.15  shows  a block  diagram  for  a phasing 
SSB  generator.  The  output  of  modulator  A is 

v (t)  = kV  cos  (u)  t)  V cos  ( a)  t).  (2.40) 

ao  m m c c 

Modulator  B is  assumed  to  be  identical  to  A,  and  its  output  is  given  by 

v (t)  = kV  cos  (w  t - tt/2)  V cos  (wt  - tt/2)  (2.41) 

bo  m m c c 


The  output  of  the  summer  is  the  lower  sideband. 


v ft)  = v (t) 
Slsb  ao 


kV  V 
c m 


COS  [ (u)  - U)  )t]  . 

c m 


(2.42) 


If  one  of  the  90°  phase  shift  networks  is  moved  to  modulator  A,  then  one  obtains  the  upper 
sideband. 


v (t)  = kV  V sin  [(u>  + u>  )t], 
usb  cm  c m 


(2.43) 


The  above  results,  however,  are  purely  theoretical.  In  practice,  one  never  achieves  perfect 
sideband  suppression  due  to  slight  variations  in  the  various  components  used.  The  most  difficult 
parameters  to  maintain  are  the  90°  phase  shift  networks.  Pappenfus  [2.27]  has  derived  formulas  for 
the  sideband  suppression  as  a function  of  the  amplitude  and  phase  imbalance.  0.5°  imbalance  in 
either  the  carrier  or  modulation  will  result  in  47  dB  sideband  suppression.  2°  gives  35  dB  and  5° 
gives  27  dB.  Likewise,  0.1  dB  amplitude  imbalance  results  in  46  dB  sideband  suppression.  0.5  dB  gives 
31  dB  and  1.0  dB  gives  25  dB . Imperfect  sideband  suppression  directly  introduces  errors  into  any 
heterodyne  measurements.  They  become  particularly  serious  when  measuring  a network  with  a sharp 
cutoff,  such  as  a bandpass  filter. 

To  receive  the  SSB  signal  and  downconvert  it  to  baseband,  a synchronous  detector  should  be  used. 

If  an  image  rejection  mixer  is  used,  additional  suppression  of  the  unwanted  sideband  is  obtained.  The 
image  rejection  mixer  is  simply  another  phasing  SSB  generator  operated  in  reverse. 

Figure  2.16  shows  a proposed  SSB  TDANA.  For  the  purposes  of  broadband  TDANA  measurements,  it  is 
desirable  to  use  a baseband  frequency  range  that  is  as  wide  as  possible.  There  are  several  conflicting 
requirements,  however.  The  major  one  is  the  bandwidth  of  the  90°  phase  shift  networks.  Most 
quadrature  (90°)  hybrids  typically  have  a bandwidth  of  one  octave.  There  are  some  broadband 
quadrature  hybrids  that  span  several  octaves,  but  their  phase  and  amplitude  specifications  are  rather 
loose.  Another  consideration  is  the  i-f  bandwidth  of  the  double  balanced  mixers  (DBM).  Most  DBM's 
typically  have  i-f  bandwidths  of  dc  up  to  1 GHz.  Some  are  available  with  i-f  bandwidths  up  to  4 GHz. 

In  light  of  these  factors  and  the  sampling  oscilloscope  and  pulse  generator  limitations,  it  appears 
that  an  i-f  frequency  of  1 to  4 GHz  would  be  optimum.  Baseband  frequencies  less  than  1 GHz  will  have 
to  be  sacrificed  in  this  arrangement  due  to  the  poor  low  frequency  performance  of  the  quadrature 
hybrids.  With  this  setup,  a single  TDANA  measurement  would  yield  data  over  a 3 GHz  frequency  span. 
Again,  if  a 10  ns  time  window  is  used,  the  frequency  spacing  of  the  resultant  data  would  be  100  MHz. 
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Figures  2.14  and  2.16  are  two  proposed  heterodyne  TDANA's.  It  is  not  clear  at  this  time  which 
system  will  give  superior  results.  Additional  study  will  be  required  before  a conclusion  can  be 
reached.  It  is  also  possible  that  a hybrid  system  may  offer  additional  advantages.  The  hybrid  would 
consist  of  an  SSB  generator  and  the  synchronous  I and  Q detectors  with  minicomputer  Hilbert  trans- 
formation. 


2.4.3  Phase  Angle  Modulation 

So  far  in  this  section  we  have  only  discussed  various  forms  of  amplitude  modulation.  The  other 
important  modulation  technique  is  angle  modulation,  either  frequency  modulation  (FM)  or  phase 
modulation  (PM).  There  is  essentially  no  difference  between  FM  and  PM  because  the  instantaneous 
frequency  oj-j  is  the  derivative  of  the  phase. 


w. 

l 


d0 

dt 


(2.44) 


or  conversely,  the  phase  is  the  integral  of  the  frequency 


{S  = / to.  dt 
i 


(2.45) 


Phase  modulation  can  be  represented  as  shown  in  figure  2.17  for  the  elementary  case  of  single-tone 
modul ation. 


v (t)  = V cos[to  t + 0(t)]  (2.46) 

pm  c c m 

0 (t)  = 0 , sin  (to  t) . (2.47) 

m d m 


0, j is  the  modulation  index.  For  this  simple  case 


v (t)  = V cos  [oo  t + 0 , sin  (to  t)] 
pm  c c d m 


(2.48) 


It  is  easily  shown  [2.25,  also  2.24,  2.28,  and  2.29]  that  Vpm(t)  is  periodic  and  may  be  expressed 
in  a Fourier  series  as 


v (t)  = V e J (0.)  cos  [(to  + nw  )t].  (2.49) 

pm  end  cm 

n=-°o 


This  function  Jn(0(j)  is  the  nth  order  Bessel  function  of  the  first  kind. 

Thus  in  sharp  constrast  to  the  AM  case  where  single  tone  modulation  created  two  sideband 
frequencies,  PM  (or  FM)  creates  an  infinite  number  of  sideband  frequencies.  In  actual  practice,  the 
Bessel  functions  are  damped  out  for  large  arguments  so  that  there  are  only  a finite  number  of 
significant  sidebands.  Nevertheless,  if  we  attempt  to  measure  an  unknown  network  or  material  with  a PM 
or  FM  test  signal,  we  end  up  receiving  a response  signal  with  many  more  unknowns  than  we  have 
measurable  quantities. 

The  effects  of  the  transmission  of  a PM  or  FM  signal  through  a linear  network  have  been 
mathematically  analyzed  [2.24].  Complete  but  extremely  involved  solutions  can  be  obtained  which 
require  computer  calculation.  Approximations  can  be  made  using  the  "quasistationary"  method  [2.24]. 
Extreme  care  must  be  exercised  in  the  use  of  this  method  because  divergent  series  are  the  typical 
output.  It  is  most  useful  when  the  transmission  characteristics  vary  in  a simple  manner. 
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Brown  has  used  PM  in  a TDANA  to  measure  the  character! sties  of  millimeter  waveguides  [2.30]. 
However,  it  is  the  opinion  of  this  author  that  PM  and  FM  are  not  valid  modulation  techniques  for 
performing  heterodyned  TDANA  measurements.  This  is  due  to  the  infinite  number  of  sideband  frequencies 
associated  with  PM  and  FM.  This  does  not  mean  that  PM  and  FM  are  of  no  value.  They  are  extremely 
useful  modulation  techniques  for  communication  purposes. 

2.4.4  Recommendations 

It  is  the  recommendation  of  this  author  that  the  most  promising  technique  for  TDANA  measurements 
in  the  millimeter  wave  region  is  the  pulse  heterodyne  method.  The  differences  between  DSB  with  I and  Q 
synchronous  detectors  and  SSB  with  an  image  rejection  mixer  need  further  detailed  study. 
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Figure  2.1  TDANA  configuration  for  S2^  measurement.  For  the  relative  positions  of  the  network 

under  test  and  the  sampler  are  interchanged. 


Voltage  (uncalibrated) 


Figure  2.2. 


Riad's  computed  impulse  response  of  the  "A-A"  Sampl 


Attenuation  in  dB 


Figure  2.3.  Riad’s  computed  frequency  response  of  the  "A-A"  Sampler. 
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Figure  2.4.  45  GHz  sine  wave  observed  on  an  "A-B"  sampler. 
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Figure  2.5.  Waveform  and  spectrum  of  an  ideal  RF  burst. 
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Figure  2.6.  Impulse  response  of  TWTAs.  Vertical:  10  V/div.  except 

(a)  which  is  3.2  V/div.  Horizontal:  200  ps  div.  (a)  NBS 

SRD  impulse  generator  (SN  3-75-3)  used  as  input  to  TWTA. 

(b)  C band,  30  dB,  10  W TWTA.  Input  impulse  attenuated 
15  dB.  (c)  X band,  30  dB,  20  W TWTA.  Input  impulse 
attenuated  16  dB.  (d)  Ku  band,  30  dB,  20  W TWTA. 
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Figure  2.7.  Spectrum  amplitudes  of  TWTA  impulse  responses. 
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Figure  2.8.  SHF  impulse  generator  pulse  output  waveform. 
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Figure  2.9.  SHF  impulse  generator  rf  pulse  output  spectrum. 
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Figure  2.11.  Simple  diode  envelope  detector. 
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Figure  2.12. 


Synchronous  detector  used  with  DSB  generator. 
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Figure  2.13.  I & Q synchronous  detector  used  with  DSB  generator. 
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Figure  2.14  Proposed  DSB  - I & Q Time  Domain 
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.15.  Phasing  type  SSB  generator. 
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Figure  2.16.  Proposed  SSB  - Time  Domain  ANA  (shown  for  USB). 


Figure  2.17.  Phase- modulation. 
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3.  IMPROVEMENT  OF  PICOSECOND-DOMAIN  SAMPLING  HEADS 


N.  S.  Nahman 
3.1  Introduction 

For  time  domain  measurements  of  repetitive  low-level  millivolt  signals  with  picosecond-domain 
resolution  (tens  of  gigahertz  bandwidth),  sampling  oscilloscopes  are  used.  The  temporal  resolution  is 
fundamentally  established  in  the  sampling  gate  structure  (i.e.,  sampling  head)  and  depends  upon  both 
the  network  and  the  gating  properties  of  the  components  comprising  the  sampling  head.  Also,  full 
utilization  of  the  inherent  sampling-head  step  response  or  temporal  resolution  requires  low  noise 
(jitter)  electronic  circuits  which  control  the  time  positions  of  the  sampling  sequence  because  time 
base  noise  introduces  a bandwidth  limitation  to  the  sampling  process  [3.1].  Automatic  pulse 
measurement  systems  can  minimize  the  effects  of  time  base  noise  and  also  provide  signal  processing  to 
enhance  the  signal -to-noi se  ratio  [3. 2, 3. 3]. 

However,  presently  available  sampling  heads  possess  several  limitations;  for  example,  their 
dynamic  range  is  limited  to  signal  levels  of  + 1 V.  Also,  at  the  instant  of  sampling,  a spurious 
strobe-ki ckout  signal  (of  several  millivolts  in  amplitude  and  several  nanoseconds  in  duration)  emanates 
from  the  signal  input  port  due  to  slight  sampling  bridge  (gate)  imbalances. 

Recently,  a new  concept  has  been  introduced  which  has  promise  to  provide  improvements  in 
sampling-gate  technology  [3.4].  The  new  concept  is  an  optically  strobed  (switched)  photoconducti ve 
sampling  gate.  This  type  of  gate  eliminates  the  problems  of  the  + 1 V dynamic  range  and 
strobe-ki ckout . 

In  p-e  time  domain  spectroscopy,  the  excitation  or  probing  signal  must  possess  a spectral 
amplitude  whose  strengths  are  great  enough  to  provide  response  signal  spectral  components  in  the  tens 
of  gigahertz  range.  There  are  two  interconnected  ways  to  visualize  the  excitation  signals.  They  are 
interconnected  because,  in  practice,  features  of  both  views  are  necessarily  present  (fig.  3.1).  On  the 
one  hand,  a suitable  probing  signal  would  be  an  impulse  of  suitable  strength,  i.e.,  one  whose  integral 
or  area  is  large  enough  (fig.  3.1a)  [3.5].  Here  the  signal  possesses  a very  short  duration 
(picoseconds)  to  insure  that  the  required  giaghertz  components  are  present.  On  the  other  hand,  a 
suitable  probing  signal  could  be  a pulse  of  non-zero  duration  and  a suitable  peak  value  (fig.  3.1b). 

For  example,  an  impulse  of  strength  S would  provide  a constant  spectrum  amplitude  of  value  S 
(fig.  3.1a).  If  the  highest  spectral  component  of  interest  must  have  at  least  a strength  of  S,  then  a 
rectangular  pulse  of  area  Ax  > S could  be  used  (fig.  3.1b).  Because  the  theoretical  impulse  function 
cannot  be  physically  realized,  non-zero  duration,  finite-magnitude  pulses  are  always  encountered  in 
practice.  Consequently,  if  the  output  magnitude  of  a given  pulse  generator  could  be  increased  by  a 
factor  of  10,  the  spectral  amplitude  would  also  increase  by  the  same  factor  (20  dB). 

With  regard  to  the  dynamic  range  of  presently  available  sampling  heads,  state-of-the-art  pulse 
generators  are  available  which  do  exceed  the  + 1 V dynamic  range.  For  example,  the  impulsive  generator 
section  of  the  NBS  SHF  pulse  generator  (delivered  to  the  USAF  Avionics  Laboratory  for  p-e  TD 
spectroscopy  appl ications)  has  a built-in  10  dB  attenuator  to  reduce  the  magnitude  of  the  impulsive 
output  so  that  it  does  not  exceed  the  sampling  head  dynamic  range.  Futhermore,  the  dynamic  range  of 
such  sampling  heads  undergoes  compression  if  the  display  is  offset  by  1 V to  take  advantage  of  the  2 V 
range  provided  by  the  + 1 V specification.  Typical  observed  compression  characteristics  for  a 10  GHz 
sinewave  are  shown  in  figure  3.2  [3.6]. 

3.2  Principles  of  Operation  for  Sampling  (Oscilloscope)  Heads 

The  sampling  gate  of  a sampling  oscilloscope  can  be  represented  by  the  equivalent  circuit  of 
figure  3.3.  The  unknown  signal  to  be  measured  is  represented  by  v-jn(t)  where  t represents  real 
time.  The  sampling  gate  is  the  variable  resistor  represented  by  the  conductance  function,  g(t).  A 
strobe,  s(t,x  ) from  the  oscilloscope  time  base  controls  the  sampling  gate  conductance.  The  instant  of 
sampling  is  represented  by  the  variable  x.  For  the  time  prior  to  x,  the  gate  is  closed,  i.e.,  g = 0. 

At  t = x,  the  strobe  opens  the  gate  allowing  current,  is(t,x),  to  flow.  This  current  deposits  charge 
on  the  sampling  capacitor,  Cs.  When  the  gate  is  again  closed,  a dc  voltage,  Vs,  remains  on  C$ 
that  is  proportional  to  an  average  of  the  signal  input  over  the  duration  of  sampling.  Vg  is  then 

measured  by  a high  impedance  voltmeter  and  displayed  on  the  cathode  ray  tube  (CRT)  screen  as  a single  dot. 
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A single  sample  is  taken  of  the  input  waveform  each  time  the  waveform  recurs.  The  time  base 
automatically  increments  t by  a small  amount  Ax  after  each  recurrence  so  that  the  waveform  is 
successively  sampled  at  a different  point  each  time  it  recurs.  The  original  waveform,  v-jn(t),  is 
thus  converted  to  an  equivalent  waveform,  Vs(x),  that  is,  stretched  into  an  equivalent  time 
variable  t. 

The  above  discussion  can  be  expressed  mathematically.  The  sampling  gate  conductance  is  assumed  to 
be  proportional  to  the  strobe  signal  and  always  of  the  same  waveshape  for  each  occurrence  (fig.  3.4). 

g(t  - x)  = ks(t  - x)  (3.1) 

where  k is  constant.  If  the  duration,  ts,  of  the  strobe  pulse  is  much  less  than  the  minimum 
circuit  time  constant,  Cs/gmax,  then  Cs  can  be  considered  to  be  a short-circuit  during  the 
sampling  interval,  ts.  The  sampling  current,  is,  is  thus  given  by 

Vt,T)  = Vin^  g^’T^'  (3.2) 


Assuming  a perfect  gate  in  the  closed  condition,  then 

i (t,x)  = 0,  for  t < x and  t > x + t . (3.3) 

s s 


The  sampling  capacitor,  Cs,  acts  as  an  integrator  such  that  the  dc  sampled  voltage,  Vs,  is 


V 

s 


1 

~ f i (t)  dt. 
C s 

s -00 


(3.4) 


Substituting  eq  (3.2)  into  eq  (3.4)  gives 


V$(x) 


~r~  / v (t)  g(t  - x)dt. 
C i n 

s -oo 


(3.5) 


For  the  special  case  when  g(t)  is  a true  impulse,  the  Vs  is  an  exact  representation  of  v-jn. 

A true  impulse  gating  function  is  the  ideal  situation  and  is  sought  after  by  sampling  gate  designers. 

In  actual  practice,  an  impulse  gating  function  can  only  be  approximated.  Equation  (3.5)  is  in  the  form 
of  the  correlation  integral;  Vs  is  the  cross  correlation  of  v-jn  and  g.  If  g(t)  is  an  even 
function,  then  eq  (3.5)  can  also  be  considered  to  be  the  convolution  integral.  For  this  analysis,  we 
will  consider  the  more  general  case  of  cross  correlation. 

As  an  example,  consider  the  case  shown  in  figure  3.4  where  s(t  - x),  and  hence  g(t  - x),  is  an 
exponential  pulse  beginning  at  t = x. 

g(t  - x)  = G exp  [-(t  - x)/o]  u(t  - x).  (3.6) 

o 


Here  u(t  - x)  represents  the  unit  step  function  beginning  at  t = x.  The  input  signal  is  assumed  to 
be  a rectangular  pulse  of  duration  tw  and  starting  at  tg; 

v.  (t)  = V [u(t  = t ) - u(t  - t - t )].  (3.7) 

in  0 o o w 


Substituting  eqs  (3.6)  and  (3.7)  into  eq  (3.5)  and  integrating  gives  the  result 
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(3.8) 


Vs(t)  = exp  [ (t  - tQ)/a],  for  T < t 


Vs'T> 


V 

s 


l-exp[ ( x-t  -t  )/a] 
U w 

l-exp(-t  /a) 
w 


• for  < T <(to  + V 


(3.9) 


V (t)  = 0,  for  t > (t  + t ) (3.10) 

s 0 w 

where 

V G a 

V$  = °c°  [1-exp  (-t w/a)]  (3.11) 

s 

Equations  (3.8)-(3.10)  are  shown  as  the  bottom  curve  of  figure  3.4. 

The  sampler  impulse  response  is  obtained  by  allowing  the  pulse  duration  tw  to  apporoach  zero  in 
which  case  eq  (3.9)  goes  to  zero.  The  impulse  response  waveform  is  thus  simply  the  mirror  image  of  the 
sampling  gate  conductance  function  waveform. 

The  sampler  step  response  is  obtained  in  a similar  manner.  In  this  case,  tw  is  allowed  to 
approach  infinity: 


Vs(t)  = Vs(step)exp[(T-tQ/a],  for  1 < tQ  (3.12) 


V$(x)  = V$(step)  = VQ Gq ct/C$ , for  y t . (3.13) 


The  sampling  efficiency  n of  the  gate  is  defined  as 


n = 


V (step) 
s 


V 


0 


(3.14) 
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3.3  Presently  Available  Devices  for  Sampling  Gates 

The  actual  sampling  gate  variable  conductance  g(t)  can  be  one  of  several  different  components 

[3.2].  At  low  frequencies  in  the  kilohertz  and  megahertz  region,  vacuum  tubes,  diodes,  and 
transistors  or  FET's  are  used.  In  the  gigahertz  and  subnanosecond  region,  commercial  sampling 
oscilloscopes  use  balanced  diode  bridges  in  2,  4,  or  6 diode  configurations  [3.7].  The  diodes  are 
normally  reverse  biased.  Upon  command  of  the  time  base,  a balanced  set  of  strobe  pulse  generators 
produce  narrow  voltage  impulses  large  enough  to  forward  bias  all  of  the  diodes.  This  has  the  effect  of 
a switch  momentarily  connecting  the  signal  input  to  the  sampling  capacitor,  Cs.  The  dynamic  range  of 
the  fast  diode  sampling  gates  is  limited  by  the  low  reverse  voltage  breakdown  of  the  diodes  and  the 
forward  contact  potential.  Severe  waveform  distortion  occurs  when  the  input  signal  is  large  enough  to 
cause  a diode  to  conduct  without  being  strobed.  Electrical  strobe  kickout  is  caused  by  slight 
imbalances  in  the  strobe  generators  and  the  bridge  circuit. 

3.4  Optically  Strobed  Sampling  Gate 

An  alternate  sampling  technique  is  an  optically  strobed  sampling  gate.  The  sampling  gate,  g(t), 
in  figure  3.1  is  a photoconductor  which  is  illuminated  by  an  optical  strobe  impulse,  s(t,e).  The  time 
resolution  achievable  with  such  a sampling  device  is  determined  by  the  duration  of  the  optical  impulse 
and  the  free-carrier  lifetime  of  the  photoconductor.  Optical  pulse  durations  of  less  than  one 
picosecond  have  been  achieved  with  a mode-locked  dye  laser  [3.8],  and  free-carrier  lifetimes  as  short 
as  50  ps  or  less  have  been  demonstrated  with  GaAs  [3.9]. 

This  type  of  sampling  device  has  the  advantage  that  the  strobe  impulse  energy  is  optical  and, 
therefore,  the  sampled  voltage  will  be  free  from  errors,  normally  found  in  the  electrically  operated 
sampling  gate,  which  are  associated  with  the  feedthrough  of  the  strobe  impulse  to  the  sampling 
capacitor.  Likewise,  there  is  no  kickout  of  the  strobe.  Also,  the  maximum  permissible  voltage  of  the 
waveform  to  be  measured  can  be  much  greater  than  is  the  case  with  conventional  sampling  oscilloscopes, 
the  dynamic  range  being  limited  only  by  voltage  breakdown  of  the  photoconductor.  In  addition,  the 
linear  V-I  characteristic  of  a photoconductor  eliminates  the  distortion  errors  found  in  an  ordinary 
sampling  gate  caused  by  diode  nonlinearities. 


3.5  Summary 

To  date,  it  has  been  demonstrated  [3.4]  that  it  is  feasible  to  use  an  optically  strobed 
photoconductor  to  sample  fast  transition,  repetitive,  electrical  transients.  Sampling  tests  have  been 
made  on  step-like  and  impulsive  waveforms.  These  tests  indicate  that  an  optically  strobed, 
photoconducti ve  gate  is  a promising  technique  for  the  measurement  of  fast  transition  waveforms  and  is 
not  limited  to  the  measurement  of  low  voltage  pulses,  as  are  present  sampling  gates.  Improvements  in 
sensitivity  and  in  transition  duration  (risetime)  can  be  expected  with  the  application  of  the 
high-power  picosecond  duration  optical  pulses  that  are  presently  available. 
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Figure  3.1a.  An  impulse,  6(t). 
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Figure  3.1b.  A rectangular  pulse. 
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Figure  3.2.  Typical  results  for  an  oscilloscope  deflection 
response  measurement  at  10  GHz. 
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Figure  3.3.  Sampling  gate  equivalent  circuit. 
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Figure  3.4.  Sampling  gate  pulse  response  waveforms. 
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4.  CORRELATION  MEASUREMENT  SYSTEMS  USING  NOISE  SIGNAL  SOURCES 


R.  A.  Lawton  and  N.  S.  Nahman 
4.1  Introduction 

The  time  domain  automatic  network  analyzer  discussed  in  Chapter  2 employs  a computer-controlled 
sampling  oscilloscope  to  acquire  the  signal  pulse  waveforms  which  are,  for  the  most  part,  deterministic 
(not  stochastic)  signals.  The  sampling  process  (equivalent  time  sampling)  is  mathematically  described 
by  the  time  integral  of  the  product  between  the  sampling  pulse  (train)  and  the  signal  [4.1].  In  this 
chapter,  a correlation  measurement  system  will  be  discussed  with  the  excitation  and  response  signals 
being  stochastic  signals.  The  excitation  source  is  a noise  generator;  consequently,  it  has  a 
randomly  varying  waveform  which  can  be  characterized  in  terms  of  statistical  measures  and  correlation 
properties. 


4.2  Correlation  Functions 

Sampling  and  correlation  measurements  have  the  same  mathematical  description,  but  they  are 
basically  implemented  by  entirely  different  methods  [4.2].  The  sampling  method  uses  a switching  or 
gating  network,!  while  correlation  in  the  picosecond  region  of  the  time  domain  employs  a nonswitched 
mixing  process  in  a nonlinear  element  to  provide  the  multiplication  process.  Both  techniques  are 
described  by  the  following  integral  equation. 


oo 

^12(t ) = f f1(t)f2(t  + T )dt 

“oo 


(4.1) 


where  x may  be  varied  continuously  or  discretely  and  in  which  t is  the  time  difference  between 
f^(t)  and  f2(t).2  Equation  (4.1)  is  commonly  called  a cross-correl ation  function  of  the  two 
functions  fq(t)  and  f2(t)  [4.3].  In  sampling,  the  sampling  pulse  or  gating  function  is 
cross-correlated  with  the  signal  waveform,  i.e.,  the  sampling  function  is  given  by 


oo 

s(t)  = / f (t)g(t  - t )dt 

-oo 


(4.2) 


where  x is  the  time  delay  of  the  gating  function  with  respect  to  the  signal  fq(t).  Here  x is 
positive  and  represents  the  sampling  point  referenced  to  the  beginning  of  the  signal  f(t). 

The  cross-correlation  function  ^ 2 1 ^ T ^ defined  by 


OO 

e>2i(x)  = / fx(t  + x)f2(t)dt. 

“00 


(4.3) 


iRefer  to  Chapter  3,  section  3.2,  for  a review  of  the  sampling  method. 


^When  both  f ]_ ( t ) and  f2(t)  are  not  square-i ntegrable  signals,  the  form  of  eq  (4.1)  is  incorrect. 
Step  functions,  random  noise,  and  periodic  signals  are  not  square-integrable  [4.3].  Under  such 
conditions,  eq  (4.1)  is  modified  to  the  form 


<12(t)  ■ lin,— 

X-*-°o 


/ f1(t)f2(t  + X )dt  . 
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Futhermore,  the  Fourier  transform  of  the  cross-correl ation  function  eq  (4.1)  is  given  by 


2(  Jco) 


°°  - It 

/ ^12(x)e  dx 


(4.5) 


= F*(ju)F2(j<o)  = F1(-ju)F2(ju) 


(4.6) 


where  F^ ( joj)  and  Fo ( juo)  denote  the  Fourier  transforms  of  f ]_ (t ) and  f2(t),  respecti vely .3 
Simi 1 arly , 


$21(jw)  = F (j(u)F  (-jto)  ; 


(4.7) 


consequently, 


$12( jw) 


(j(*>). 


(4.8) 


4>12(ju))  and  $2i(jco)  are  complex;  each  is  referred  to  as  a cross-correl  ation 

spectral-density  function.  When  fi(t)  and  f2(t)  are  identical,  eq  (4.1)  yields  the  autocorrelation 
function  of  fj (t) 


VT) 


00 

J f1(t)f1(t  + x )dt . 

_co 


(4.9) 


The  Fourier  transform  of  the  autocorrelation  function  eq  (4.9)  yields  | F^(ju))|  2 which  is  a real 
and  positive  function  of  co.  | F^ ( jco )|  2 is  commonly  called  the  energy-spectral-density  function; 
it  is  related  to  the  spectrum  amplitude  A(co)  through  the  relation  [4.4] 

A(u)  « 2 |F1(ju)|.  (4.10) 


■3 1 f 0i2(x)  is  periodic,  the  integral  form  must  be  modified  as  discussed  in  footnote  2. 


4.3  Using  Noise  and  Correlation  Methods  for  p-e  Measurements 


An  alternate  technique  for  determining  network  characteristics  and  thus  p and  e over  a wide  band 
of  frequencies  is  the  noise  cross-correlation  approach.  The  basic  idea  of  this  approach  can  be 
introduced  with  the  autocorrel ator  outlined  in  figure  4.1  which  could  be  implemented  in  any  desired 
waveguide,  e.g.,  coaxial  line,  WR12  (60-90  GHz),  etc.  As  shown  by  the  diagram,  a noise  signal  with  a 
time  domain  waveform,  f(t),  is  split  into  two  channels,  one  of  the  channels  being  delayed  with  respect 
to  the  other.  The  two  signals  are  then  recombined,  and  the  result  multiplied  and  integrated  which 
results  in  the  function  A + 20n(x)  at  the  output  of  the  integrator, 


A + 2<(>^  (t  )- 


lim  2T 
T 


/ 


fl(t) 


-T 


f (t+  r)dt . 


(4.11) 


The  square  law  device  produces  the  output 

g(t)  = f*(t  - Tx)  + f^(t  - x)  + 2 f (t  - T^f^t  - t) 


which  the  integrator  transforms  to 

/g(t)dt  = < f*(t)  > + < f*(t)  > + 2 / f1(t)f1(t  + T)dt  . 


(4.12) 


(4.13) 


2 

The  first  two  terms  are  identical  and  are  each  equal  to  the  time  average  value  of  fj(t);  the  third 
term  is  twice  the  autocorrelation  function,  0n(x),  of  fi(t).  Consequently,  the  output  of  the 

integrator  is  20n(x)  plus  a constant,  A.  The  signal  path  delays,  t]_  and  x,  are  such  that  the 

two  signals  can  be  in  effect  delayed  or  advanced  with  respect  to  each  other. 

An  example  of  an  experimentally  measured  autocorrelation  function  for  a picosecond  domain  pulse 
signal  is  shown  in  figure  4.2  [4.3].  Because  the  autocorrelation  function  is  always  an  even  function 
of  the  time  delay,  only  the  portion  for  f 0 is  shown. 

Let  us  now  consider  figure  4.3  which  is  exactly  the  same  as  in  figure  4.1  except  that  a network 
whose  impulse  response,  h(x),  is  to  be  determined  has  been  inserted  in  the  upper  channel.  Now  at  the 

output  of  the  integrator,  apart  from  a constant,  A2,  one  has  the  cross  correlation,  2021 (t), 

between  fi(t)  and  the  output  of  the  unknown  network. 

1 T 

4,  (x)  = lim  2T  / f2(t)  f x (t  + T)dt.  (4.14) 

-T 

Here  the  subscripts  1 and  2 denote  the  input  and  output  functions  of  the  unknown  network.  Note  that 
f2(t)  is  the  convolution  between  f 1 ( t ) and  h(t)  denoted  by  fi(t)*h(t),  i.e., 

t 

f2(t)  = fj  (t)  * h(t)  = f f^t  - X)  h(A)dX  . (4.15) 

0 


The  autocorrelation  function  will  now  be  denoted  as  0n(x). 

In  reference  [4.4]  it  is  seen  that  the  Fourier  transform  of  0n(x)  yields  the  input  power 
density  spectrum,  Su(u)),  of  the  signal,  f(t);  and  the  Fourier  transform  of  02l(x)  yields 
the  cross  spectrum,  S2 1 ( co) . The  significance  of  the  cross  spectrum  is  that  when  divided  by  the 
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input  power  density  spectrum,  Sn(w),  it  yields  the  complete  system  transfer  function,  H(u>),  of  the 
inserted  network.  In  other  words, 


S21M/Sn(u))  = H(w). 


(4.16) 


Now  let  us  consider  a system  in  which  the  signal  source  is  a white  noise  generator  whose  power 
density  spectrum  has  the  value  K (independent  of  frequency).  In  this  case 

S21(oj)  = KH(w)  (4.17) 


which  has  the  inverse  transform 

*21(t)  = 2mKh(T).  (4.18) 


In  other  words,  for  a white  noise  source,  the  system  impulse  response  is  directly  proportional  to 
the  cross-correl ati on  function  between  the  input  and  output  signals.  From  this  impulse  response  or, 
equivalently,  from  the  system  transfer  function,  all  the  network  parameters  of  the  inserted  network 
such  as  dielectric  constant  and  permeability  could  be  determined. 

While  the  above  description  has  been  given  in  general  for  any  waveguide,  it  can  be  applied  to  any 
specific  waveguide  band.  However,  every  band  will  be  somewhat  band  limited,  and  the  noise  source,  for 
that  reason  alone,  will  not  be  a completely  white  noise  source.  Equation  (4.16)  can  then  be  used  to 
remove  the  effect  of  the  band-limited  source.  Some  examples  of  white  noise,  cross-correlation 
functions  or  network  impulse  responses  that  one  might  encounter  are  those  of  the  low  pass  R-C 
integrator  and  the  Debye  dielectric  lossy  transmission  line  [4.5]  which  are  given  in  figure  4.4  along 
with  their  respective  equivalent  circuits. 

4.3.1  Time  Resolution  Requirements 

The  next  question  one  might  ask  is  what  time  delay  temporal  resolution  is  required.  To  get  a 
rough  estimate,  suppose  that  the  filtering  characteristic  of  the  inserted  network  is  approximately 
Gaussian  in  character;  then  the  transfer  function  H(jw)  would  be  approximated  by 

H(jw)  = e"“  /4k  , (4.19) 


while  the  corresponding  approximate  network  impulse  response  would  be  the  inverse  Fourier  transform 
of  H(jto)  which  is 

2 2 

h(  t)  = ~~  e"k  (4.20) 

2tt 

The  two  functions  H(jto)  and  h(x)  are  each  Gaussian  functions  and  possess  the  1/e  intervals  of 

Aw  = 4k  (4.21) 

and 

At  = 2/k , (4.22) 

respectively  (fig.  4.5).  Therefore,  for  a network  bandwidth  of 

Aw 

Af  = — — = 30  GHz  , 

2tt 
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which  would  totally  occupy  the  available  bandwidth  in  the  UR-12  waveguide  band  (60-90  GHz),  the 
corresponding  impulse  response  duration  would  be 


At 


2 

1 5ttx10  9 


21.2  ps. 


This,  in  turn,  would  require  that  the  delay  measurement  resolution  be  of  the  order  of  two  picoseconds 
or  better,  which  corresponds  to  a distance  measurement  of  0.6  mm  for  electromagnetic  wave  propagation 
in  air. 


4.3.2  Other  Realization  Considerations 

Square  law  detector:  Diodes,  thermistors,  and  resistive  bolometers  with  pyroelectric  sensors  may 

be  used  to  provide  the  multiplying  function.  The  use  of  resistive  heaters  and  pyroelectric  sensors  is 
outlined  in  detail  in  [4.3].  The  same  reference  describes  how  the  device  capacitance  can  be  used  for 
the  required  integration.  The  use  of  diodes  can  provide  greater  sensitivity  than  the  use  of  bolometers 
or  pyroelectric  detectors,  but  the  latter  two  devices  provide  better  square  law  response  with  the 
pyroelectric  detector  being  the  more  sensitive  of  the  two. 

Noise  source:  The  solid  state  noise  source  developed  by  Kanda  [4.6]  is  a strong  source  and  thus 

provides  large  noise  temperatures , typically  290,000  K. 
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Figure  4.2.  Autocorrelation  function  as  measured. 
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Figure  4.4a.  Impulse  response  of  low-pass  RC  network.  Figure  4.4b.  Impulse  response  of  a Debye 

dielectric  lossy  transmission 


Figure  4.5.  Gaussian  system  responses. 
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5.  THE  FREQUENCY  DOMAIN  SIX-PORT  AUTOMATIC  NETWORK  ANALYZER 


C.  A.  Hoer  and  C.  M.  Allred 
5.1  Introduction 

The  six-port  concept  has  been  successfully  applied  in  the  realization  of  microwave  reflectometers 
[5. 1,5. 2]  and  vector  voltmeters  [5.3].  A conventional  six-port  refl ectometer  is  a six-port  junction 
designed  to  measure  reflection  coefficient  and  power  at  one  of  its  ports  in  terms  of  power  readings 
taken  at  four  of  the  other  ports  when  an  rf  signal  is  applied  at  the  remaining  port.  The  measured  or 
observed  quantities  are  scalar  as  opposed  to  phasor  quantities,  i.e.,  the  observables  are  ratios  of 
transmitted  power  at  the  ports  rather  than  complex  voltages  or  currents.  A typical  realization  of  a 
six-port  refl ectometer  is  shown  in  figure  5.1  [5.4]. 

The  value  of  the  complex  reflection  coefficient  p obtained  from  a six-port  refl ectometer  is 
calculated  from 


P 


Vc.P.  + jVs.P. 
Va.P. 


(5.1) 


which,  with  few  exceptions,  applies  to  any  six-port  junction  [5.2].  The  P^'s  are  the  power  readings 
from  four  of  the  five  detectors  in  figure  5.1a,  where  P3  must  be  one  of  the  detectors  used.  The  12 
real  constants,  c-j , s-j , and  a-j , characterize  the  six  port.  They  are  functions  only  of  the 
scattering  parameters  of  the  six  port  and  of  the  reflection  coefficient  of  the  detectors.  They  are 
independent  of  the  properties  of  the  rest  of  the  measurement  system  including  the  flexible  cables  and 
the  generator.  One  can  divide  numerator  and  denominator  of  eq  (5.1)  by  any  one  of  these  12  constants 
leaving  11  new  constants  which  must  be  determined  by  some  calibration  process  such  as  those  mentioned 
in  section  5.3. 


Using  all  five  detectors,  one  can  display  approximate  values  of  p on  an  oscilloscope  as  shown  in 
figure  5.1b.  This  application  is  based  on  the  following  equation  for  this  particular  seven-port 
junction  assuming  ideal  (i.e.,  lossless,  infinite  isolation,  zero  reflection,  etc.)  components  [5.2]: 


P 


(P  - P ) + j(P  - P ) 
v 5 l'  6 8 


(5.2) 


As  this  equation  shows,  P3  acts  like  a scale  factor.  P3  can  be  measured  to  determine  the  scale 
of  | p | on  the  oscilloscope  screen.  It  can  also  be  used  to  control  the  level  of  the  generator.  A 
real-time  analog  display  of  p is  especially  useful  in  checking  out  the  measurement  system. 

By  using  two  six-ports  a network  analyzer  can  be  built  which  can  be  calibrated  using  a length  of 
precision  transmission  line  as  the  only  standard  needed.  Two  six-port  reflectometers  are  used  in  the 
configuration  shown  in  figure  5.2  to  measure  all  the  network  parameters  of  any  linear  reciprocal  or 
nonreciprocal , active  or  passive  two-port,  inserted  between  the  two  reflectometers.  This  configuration 
is  most  similar  to  the  interference  bridge  described  in  1962  by  Altschuler  [5.5].  The  two  measurement 
systems  differ  in  that  six-port  reflectometers  are  used  in  place  of  slotted  lines,  isolators,  and 
tuners.  This  change  makes  the  present  system  easy  to  automate  and  broadband,  in  the  sense  that 
measurements  can  be  made  at  any  fixed  frequency  from  2 to  18  GHz  after  proper  calibration  at  that 
frequency.  Four-port  reflectometers  could  be  used  in  place  of  six-port  ref  1 ectometers , but  the  system 

would  not  be  as  simple  as  the  one  described  here.  A great  advantage  of  the  six-port  is  that 
imperfect,  and  hence  very  broad  band  components,  Q and  D in  figure  5.1,  are  used;  all  imper- 
fections are  accounted  for  in  the  calibration  rather  than  tuned  out  at  successive  frequencies. 

The  measurement  system  described  here  differs  from  most  other  network  analyzers  in  that  the 
six-port  reflectometers  can  be  made  small  enough  and  portable  enough  to  connect  directly  to  the  device 
under  test,  eliminating  flexible  cables  or  arms  between  the  measurement  system  and  the  device  being 
measured.  Flexible  cables  "behind"  each  six  port  do  not  enter  into  the  system  calibration  except  in 
the  measurement  of  the  phase  shift  through  nonreci procal  two  ports.  The  present  system  is  simpler  than 
most  other  netowrk  analyzers  in  that  no  i-f  source  is  used  and  no  phase  detectors  are  required. 
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Self-calibration  techniques  similar  to  those  described  by  Allred  and  Manney  [5.6]  for  calibrating 
two  four-port  couplers  can  be  used  to  calibrate  two  six-port  ref 1 ectometers . When  applied  to  the  dual 
six-port  measurement  system,  this  technique  yields  the  constants  in  eq  (5.1)  needed  to  calculate  pj 
and  P2»  Only  one  standard  is  required  in  the  self-calibration  technique.  This  standard  can  be  a 
uniform  length  of  transmission  line  or  one  known  termination.  Details  of  this  self-calibration 
technique  are  given  in  references  [5. 7, 5. 8]. 


5.2  Theory  of  the  Dual  Six-Port  Network  Analyzer 
A.  Determination  of  S22>  | S^2  |>  | S21  | 

The  basic  setup  for  using  two  six-port  ref 1 ectometers  to  measure  all  of  the  scattering  parameters 
of  a two  port  is  shown  in  figure  5.2.  The  rf  signal  is  applied  to  both  ref 1 ectometers  simultaneously. 
Ref 1 ectometer  1 measures  the  complex  ratio  pj  = b^/a2  at  reference  plane  1.  Reflectometer  2 
measures  P2  = b2/a2  at  reference  plane  2.  pi  and  and  p2  are  not  reflection  coefficients 

in  the  normal  sense,  but  simply  relate  two  traveling  waves  in  opposite  directions.  These  two  ratios 
are  related  by  the  S-parameters  of  the  two  port  inserted  betwen  the  two  ref 1 ectometers.  Let  the 
scattering  equations  for  the  two  port  under  test  be 


bl  = 


S a 
11  1 


+ S,  „a„ 
12  2 


(5.3) 


b = S a + S a 
2 21  1 22  2 


(5.4) 


Dividing  eq  (5.3)  by  aj  and  eq  (5.4)  by  a2  gives 

b 


_ _J_  _ , __JL 

P1  ax  S11  S12  a1 


P2  " 77  “ S22  * S21  77 

2 2 


(5.5) 


(5.6) 


Since  | a2/ai  | can  be  greater  or  less  than  one,  | p]_  | and  | P2  | can  also  be  greater  or  less  than 

one.  Six-port  ref 1 ectometers  are  capable  of  measuring  values  of  | p | both  larger  and  smaller  than  one. 


Eliminating  a2/a^  from  eqs  (5.5)  and  (5.6)  gives 


(pi  - sn)  (p2  - s22)  - s12s21 


(5.7) 


or 


32S11 


+ P1S22  ' A = 


P1P2 


(5.8) 


where  A is  the  determinant  of  the  scattering  matrix. 


A 


= S S 
11  22 


s s 

12  21 


(5.9) 


Three  equations  like  eq  (5.8)  are  solved  for  S^,  S22,  and  A.  These  equations  are  generated 
by  measuring  p^  and  P2  for  three  different  values  of  a2/aq  which  are  determined  by  the  setting 
of  the  attenuators  l\\  and  A2  and  of  the  phase  shifter  0.  The  values  of  A} , A2,  and  0 do 
not  need  to  be  known  since  they  do  not  enter  into  the  equations. 
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If  it  is  known  that  the  two  port  under  test  is  reciprocal,  so  that  = S21,  the 
magnitude  can  be  obtained  from  eq  (5.9), 


or  from  eq  (5.7) , 


21 


= lsns22 


(pi  - snHp2  - s22)  |. 


(5.10) 


(5.11) 


After  three  or  more  equations  like  eq  (5.8)  ar’e  solved  for  Sqq,  S22,  and  A,  then  Si 1 and 
S22  can  be  used  in  eq  (5.11)  to  calculate  values  of  | Sq2  | and  | S2q  | for  each  measurement. 

The  scatter  in  the  values  of  | S^2  | gives  an  indication  of  the  precision  with  which  I S^2  | has 
been  measured. 

If  it  is  not  known  that  the  two  port  under  test  is  reciprocal,  | S}2  | and  | S2q  | are 
obtained  from  eqs  (5.5)  and  (5.6)  which  give 


12 


(5.12) 


(5.13) 


The  six-port  ref 1 ectometers  can  be  used  to  measure  | a2/aq  | as  well  as  pq  and  P2  in  these 
equations.  Altschuler  [5.4]  discusses  the  advantages  of  making  | a2/aq  | large  or  small  to  increase 
the  sensitivity  when  measuring  small  values  of  Sq2  or  S2q. 


B.  Determination  of  i!)q2  = ^21 


The  equations  given  in  sections  5.2A  can  be  used  to  determine  all  of  the  scattering  parameters 
except  ip 12  and  ij^q,  the  phase  angles  of  Si 2 and  S21 • For  either  a reciprocal  or 
nonreciprocal  two  port,  eqs  (5.5)  and  (5.6)  give 


K12 


*1  ' 


(5.14) 


21 


'Pi  + 'P. 


(5.15) 


where 


^ 5 arg(Pl  - Sn) 

(5.16) 

\P2  = arg(p2  - S22) 

(5.17) 

ip  = arg(a  /a  ). 

6 Lx. 

(5.18) 

Since  the  six-port  ref 1 ectometers  do  not  directly  measure  the  phase  angle  of  a2/ai,  the  angle 
'Pa  is  not  yet  known.  The  angles  ^q  and  1P2  are  known  from  measured  values  of  pq  and  p2  and 
computed  values  of  Sqq  and  S22  (see  sec.  5.2A). 
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When  the  two  port  is  reciprocal,  eqs  (5.14)  and  (5.15)  can  be  solved  for  both  >la  and 
i>\2  = ^21 5 


ip  + w 

V1  v2 


^12  ^21  2 


+ nrT  , n = 0,  1 


(5.19) 


i(j1  - tp2 

W = - — + rnr  , n = 0,  1 

a l. 


(5.20) 


where  the  ambiguity  rnr  airses  from  the  division  by  2,  and  n is  an  integer.  All  phase  angles 
considered  here  are  modulo  2tt  as  in  any  frequency  domain  measurements.  It  can  be  seen  from  eqs  (5.14) 
and  (5.15)  that  n must  take  on  the  same  value  in  eqs  (5.19)  and  (5.20).  The  value  of  n can  be 
determined  from  either  an  estimate  of  ^32  or  ipa.  A good  estimate  of  \pg  can  be  obtained  from 
equations  derived  for  the  nonreciprocal  case  which  is  considered  in  section  5.2C. 

The  measurement  technique  described  up  to  this  point  is  similar  to  methods  which  use  only  one 
ref 1 ectometer  to  measure  P3  when  the  two  port  is  terminated  with  three  different  known  terminations 
at  reference  plane  2.  The  equations  which  apply  when  using  known  terminations  are  the  same  as 
eqs  (5.3)  to  (5.20)  if  p£  is  replaced  by  1/T[_  where  T[_  is  the  reflection  coefficient  of  the 
termination  at  refernce  plane  2.  The  difference  here  is  that  a2  can  be  made  much  larger  than  if 
passive  terminations  are  used.  This  makes  it  possible  to  measure  larger  values  of  attenuation  using 
the  setup  in  figure  5.2.  Another  difference  is  that  the  ratio  | a2/ai  | cannot  be  measured  when 
passive  terminations  are  used  so  that  S32  and  S21  of  nonreciprocal  two  ports  cannot  be 
measured. 

C.  Determination  of  ^32  i ^21 

When  the  two  port  under  test  is  not  reciprocal,  ^32  and  ^21  must  be  calculated  from 
eqs  (5.14)  and  (5.15)  where  j;a  can  be  determined  as  follows.  The  determination  of  depends  in 
part  on  knowledge  of  the  measurement  system  which  can  be  viewed  as  a three-port  junction  whose  ports 
are  1,  2,  and  3 (near  the  generator)  in  figure  5.2.  For  this  three-port  junction,  one  can  show  that 


where  the  small  s-jj's  are  the  scattering  parameters  of  the  equivalent  three-port  measurement 
system.  Rewrite  eq  (5.21)  as 


a 


1 


Clpl 


(5.22) 


where  C3,  C2,  and  C3  are  defined  by  corresponding  terms  in  eqs  (5.21)  and  5.22).  The  C's  for  a 
particular  setting  of  A3,  A2,  and  <t>  can  be  determined  from  three  equations  like  eq  (5.22)  when 
PI,  P2,  and  a2/ai  are  known,  which  is  the  case  when  reciprocal  two  ports  are  measured,  solving 
eq  (5.22)  for  82/83, 


C3  + C1P1 
1 + C2P2 


(5.23) 


Once  the  C's  are  known,  ij;a  can  be  calculated  from  eq  (5.23)  when  measuring  either  reciprocal  or 
nonreciprocal  two  ports. 
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When  the  two  port  is  nonreciprocal,  calculated  from  eq  (5.23)  is  used  in  eqs  (5.14  and  (5.15) 
to  calculate  i>i2  and  471  • When  the  two  Port  ’’s  reciprocal,  the  value  of  4*a  from  eq  (5.23)  is 
considered  to  be  only  an  estimate  for  determining  n in  eq  (5.20).  This  is  because  the  value  of  4;a 
calculated  from  eq  (5.20)  will,  in  general,  be  more  accurate  than  that  calculated  from  eq  (5.23). 

Since  the  C's  in  eq  (5.23)  are  functions  of  the  scattering  parameters  of  the  three-port  measurement 
system,  they  are  functions  of  the  stability  of  the  flexible  cables  and  also  of  the  setting  and 
repeatabi 1 ity  of  the  components  l\\,  A2,  and  0.  None  of  the  other  equations  up  to  eq  (5.21)  are 
functions  of  these  components  or  the  cables.  For  this  reason,  eq  (5.23)  is  used  only  to  calculate  ipa 
and  not  | a 2/ai  | which  can  be  obtained  from  the  two  six-port  reflectometers  independently  of  the 
properties  of  the  three-port  measurement  system.  Then  all  of  the  calculated  S-parameters  of  the  two 
port  under  test  will  be  independent  of  the  properties  of  the  three-port  measurement  system  except 
ipl2  and  471  of  nonreciprocal  two  ports. 

5.3  Using  Six-P-ort  Systems  for  p-e  Measurements 
5.3.1  The  Six-Port  Ref 1 ectometer 

One  way  of  measuring  the  magnetic  and  dielectric  properties  of  small  samples  at  fixed  frequencies 
is  with  a six-port  reflectometer  (fig.  5.3).  A six-port  reflectometer  is  a six-port  junction  such  as 
shown  in  figure  5.1  designed  to  measure  reflection  coefficient  and  power  at  one  of  its  ports  in  terms 
of  power  readings  taken  at  four  of  the  other  ports  when  an  rf  signal  is  applied  at  the  remaining  port. 
The  complex  reflection  coefficient,  r,  is  given  in  terms  of  the  four  power  readings,  P3...P5,  by 
eq  (5.1)  or  more  compactly  by 


E z . P . 

7 V .1=3,  4,  5,  6,  (5.24) 


where  the  four  complex  constants  z-j  and  the  four  real  constants  a-j  are  functions  of  the  six-port 
junction.  These  constants  are  obtained  by  a calibration  process  involving  at  least  one  but  usually 
more  known  terminations  which  are  connected  to  the  measurement  port.  The  properties  of  a sample  can 
be  obtained  from  measurements  of  r with  and  without  the  sample  in  the  holder. 

The  main  advantage  of  the  six-port  reflectometer  over  the  conventional  four-port  reflectometer  is 
that  the  phase  angle  as  well  as  magnitude  of  r is  obtained  from  simple  power  measurements.  It  is  not 
necessary  to  have  a phase  detector  on  the  sidearem  of  the  six  port  to  measure  the  phase  of  r.  Only  one 
source  is  used  since  no  frequency  conversion  is  involved.  The  six  port  is  relatively  easy  to  automate 
since  all  information  is  obtained  from  a set  of  dc  voltage  readings  taken  from  the  four  power  meters. 

The  resolution  in  measuring  r is  essentially  the  resolution  in  measuring  the  four  values  of  power 
which  is  one  part  in  10^  for  thermistor  types  of  power  detectors.  These  detectors  are  commercially 
available  up  to  100  GHz.  Waveguide  components  for  constructing  six-port  junctions  are  also 
commercially  available  to  100  GHz.  Therefore  it  should  be  possible  to  construct  six-port 
reflectometers  in  each  waveguide  band  up  to  100  GHz  and  possibly  even  higher. 

An  experimental  automated  waveguide  six-port  reflectometer  has  been  constructed  at  NBS  for 
measurement  of  power  and  T in  the  50  to  75  GHz  range.  Preliminary  data  indicate  that  it  has  a 
resolution  of  0.0001  in  measuring  r in  this  frequency  range.  The  accuracy  has  not  yet  been  determined, 
but  for  small  changes  in  r(A  IX  0.1),  the  accuracy  of  the  change  is  expected  to  be  limited  primarily  by 
the  resolution.  The  time  required  to  measure  r is  several  seconds  per  frequency.  Most  of  that  time  is 
required  to  set  and  phase-lock  the  frequency.  This  time  can  be  reduced  considerably  by  not  locking  the 
frequency,  but  the  accuracy  would  be  decreased. 

5.3.2  The  Dual  Six-Port  Network  Analyzer 

Another  way  of  measuring  the  properties  of  a sample  is  to  insert  it  in  a two-port  sample  holder 
which  is  inserted  between  two  six-port  reflectometers  as  shown  in  figure  5.4.  The  two  six-port 
reflectometers  are  capable  of  measuring  all  of  the  scattering  parameters  of  the  two  port  with  and 
without  the  sample  inserted  in  the  holder.  These  measurements  should  yield  both  the  transmission  and 
reflection  properties  of  the  sample. 
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There  are  several  advantages  in  having  two  six-port  reflectometers  instead  of  one.  It  is  much 
easier  to  calibrate  two  six-port  reflectometers  than  it  is  to  calibrate  just  one  in  that  fewer 
connections  need  to  be  made  by  the  operator,  and  fewer  standards  are  required.  The  only  standards 
needed  to  calibrate  a pair  of  six  ports  are  a flat  short  and  a length  of  precision  waveguide.  A 
typical  calibration  procedure  is  outlined  in  figure  5.5 

Another  advantage  is  that  by  connecting  the  two  six-port  measurement  ports  together,  one  can 
determine  if  the  two  reflectometer  measurements  of  V agree  with  each  other.  The  stability  of  the 
calibration  with  time  can  thus  be  determined  and  monitored. 

A pair  of  coaxial  six-port  ref lectometers  have  bveen  constructed  at  NBS  for  measuring  S-parameters 
of  two-port  devices  in  the  range  of  2 to  18  GHz.  The  techniques  and  software  for  calibrating  this  pair 
of  reflectometers  should  be  applicable  to  waveguide  systems  at  higher  frequencies.  The  precision  of 
the  NBS  dual  six-port  system  in  measuring  the  S-parameters  of  a two  port  is  approxi amtely  0.0001. 

One  disadvantage  in  measuring  all  of  the  S-parameters  is  that  it  takes  about  eight  times  as  long 
as  it  does  to  measure  T with  a single  six  port. 
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(b) 


A seven-port  ref lectometer  constructed  with  quadrature  hybrids  Q, 

and  in-phase  power  dividers  D,  in  a conventional  circuit  to  measure 

reflection  coefficient  P.  A six-port  ref lectometer  is  obtained  by 

simply  ignoring  one  of  the  four  outputs  P ...  P . If  the  detectors 

3 o 

each  have  an  output  dc  voltage  proportional  to  the  input  RF  power,  the 
voltages  can  be  applied  to  an  oscilloscope  as  in  (b)  to  get  a display 
of  reflection  coefficient  P . The  constant  k need  not  be  known  to 
display  T.  When  the  termination  is  replaced  with  a two  port,  the  ratio 
p is  displayed  instead  of  T. 

Figure  5.1.  A six-port  ref lectometer . 
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Figure  5.2.  Basic  circuit  for  measuring  all  of  the  scattering 
parameters  of  a two  port  with  two  six-port  ref lectometers . 
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Figure  5.3.  Six-port  ref lectometer  y-e  measurement  system. 
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Figure  5.4. 


The  dual 


six-port  network  y-e 


measurement 


system. 
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DUAL  6-PORT  ANA 


'■6-PORT  TSD"  CALIBRATION  TECHNIQUE 


Figure  5.5.  The  four  different  connections  made  in  calibrating  the  dual  6-port 
ANA  with  the  "6-port  TSD"  technique,  where  TSD  stands  for  "Thru",  "Short",  "Delay". 
The  parameters  of  the  normal  10  dB  pad  do  not  need  to  be  known. 
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6.  CONCLUSION 


6.1  Summary 

The  purpose  of  this  report  was  to  describe  new  methods  for  implementing  automatic  p-e  measurements 
of  radar  absorber  materials  for  applications  primarily  above  20  GHz.  After  reviewing  the  literature 
and  drawing  upon  the  collective  experience  of  the  NBS  staff,  three  methods  were  selected  for 
discussion: 

1.  The  Time  Domain  Automatic  Network  Analyzer. 

2.  Correlation  Measurement  Systems  Using  Noise  Signal  Sources. 

3.  The  Six-Port  Frequency  Domain  Automatic  Network  Analyzer. 

Because  there  were  several  technical  options  in  implementing  a time  domain  automatic  network 
analyzer,  consideration  had  to  be  given  to  four  topics: 

a)  Extension  of  present  baseband  system  performance  to  higher  frequencies. 

b)  Millimeter  wave  pulse  generators. 

c)  Heterodyne  techniques. 

d)  Improvement  of  picosecond-domain  sampling  heads. 

For  each  of  the  three  new  methods  discussed  in  this  report,  related  time  domain,  correlation  and 
microwave  research  and  development  have  been  implemented  at  NBS  in  the  past  for  other  applications. 
Consequently,  all  three  methods  are  feasible. 


6.2  Discussion 

The  central  question  to  be  addressed  is  simply,  "which  method  or  methods  would  achieve  the  most 
effective  results?" 

There  are  two  important  factors  to  consider  in  answering  the  question.  The  first  is  the  time  and 
effort  required  to  refine  the  method(s)  for  application  above  20  GHz  and  into  the  millimeter  wave 
region  to  about  100  GHz.  The  second  factor  depends  upon  the  sample-holder  structure. 

In  the  technical  descriptions  presented  in  Chapters  2,  4,  and  5,  it  was  assumed  that  the  sample 
fits  exactly  within  the  system  waveguide,  i.e.,  the  sample  holder  was  a uniform  piece  of  waveguide  which 
was  completely  filled  by  the  sample,  and  the  sample  faces  were  plane  and  parallel  to  the  wavefront  of 
the  principal  waveguide  mode.  These  assumptions  become  increasingly  difficult  to  satisfy  when: 

(1)  waveguide  dimensions  are  decreased  to  allow  only  principal  mode  propagation,  and  (2)  the  sample 
physical  properties  are  such  that  the  sample  is  difficult  to  shape  and  susceptible  to  damage  which  can 
alter  its  electromagnetic  and/or  physical  properites.  Consequently,  these  assumptions  become 
increasingly  untenable  with  increasing  frequency.  The  only  alternative  available  is  to  depart  from  the 
principal  mode  sample  holder  and  employ  an  electromagnetically  large  or  over-moded  sample  holder. 

Physically  large  sample  holders  could  approach  free-space  conditions  in  which  the  sample  is  placed 
between  two  antennas  and  irradiated  by  a characterized  wavefront  eminating  from  one  of  the  antennas. 

The  wavefront  passes  through  the  sample  and  is  received  by  the  second  antenna.  However,  such  a system 
would  ideally  require  a very  large  sample,  i.e.,  a large,  thin,  plane  sheet. 

In  practice,  available  samples  are  10  to  15  cm  (4  to  6 inches)  square.  Consequently,  if  such 
samples  are  to  be  directly  used  without  physical  alteration,  the  sample  holder  would  have  to  be  one 
which  would  accept  the  whole  sample  and  provide  a charcterizable  electromagnetic  coupling  to  the 
sample.  Clearly,  the  sample  holder  would  be  over-moded  and  for  the  higher  frequencies  approach  a quasi 
free-space  condition.  The  critical  factor  which  must  be  resolved  for  a successful  application  of  an 
over-moded  or  quasi  free-space  sample  holder  is  the  characterization  of  the  electromagnetic  coupling  to 
the  sample.  If  the  coupling  is  known,  then  transmission  and/or  reflection  data  can  be  analyzed  to 
yield  the  electromagnetic  parameters  p and  e. 
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It  is  important  to  recognize  that  both  the  time  domain  and  frequency  domain  p and  e measurement 
methods  require  characteri zation  of  the  over-moded  electromagnetic  coupling  to  the  sample. 
Consequently,  this  means  that  the  network  properties  or  physical  model  must  be  determined  to  establish 
the  transfer  functions  between  the  uniform  waveguide  reference  planes  and  the  sample.  Neither  the 
frequency  domain  nor  the  time  domain  approach  eliminates  the  need  to  establish  the  transfer  functions 
of  the  over-moded  or  quasi  free-space  sample  holder. 

6.3  Recommendations 

In  answer  to  the  central  question  stated  at  the  beginning  of  section  6.2,  it  is  recommended  that 
future  work  include  the  following  technical  approaches: 

1.  Heterodyne  Time  Domain  Automatic  Network  Analyzer. 

2.  Six-Port  Automatic  Network  Analyzer. 

3.  Over-moded  sample  holders. 

The  first  two  approaches  primarily  require  a development  effort  utilizing  available  microwave  or 
millimeter  wave  components  and  available  mini -computer  systems  and  interfaces.  The  third  approach 
requries  a research  effort  to  establish  the  properties  of  over-moded  or  quasi  free-space  sample 
hoi ders . 
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